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ABSTRACT 

We present results based on an XMM-Newton observation of the high luminosity narrow-line QSO 
PHL 1092 performed in 2003 January. The 0.3 — 10 keV spectrum is well described by a model which 
includes a power-law (r ~ 2.1) and two blackbody components (kT ~ 130 eV and kT ~ 50 eV). The 
soft X-ray excess emission is featureless and contributes ~ 80% to the total X-ray emission in the 
0.3 — 10 keV band. The most remarkable feature of the present observation is the detection of X-ray 
variability at very short time scale: the X-ray emission varied by 35% in about 5000 s. We find that 
this variability can be explained by assuming that only the overall normalization varied during the 
observation. There was no evidence for any short term spectral variability and the spectral shape 
was similar even during the ASCA observation carried out in 1997. Considering the high intrinsic 
luminosity (~ 2 x 10 45 erg s _1 ) and the large inferred mass of the putative black hole (~ 1.6x 10 8 M Q ), 
the observed time scale of variability indicates emission at close to Eddington luminosity arising from 
very close to the black hole. We suggest that PHL 1092 in particular (and narrow line Seyfert galaxies 
in general) is a fast rotating black hole emitting close to its Eddington luminosity and the X-ray 
emission corresponds to the high-soft state seen in Galactic black hole sources. 
Subject headings: galaxies: active — galaxies: individual (PHL 1092) — X-rays: galaxies 



1. INTRODUCTION 

Narrow Line Seyfert 1 (NLS1) ga laxies were first clas- 
sified as a subclass of Seyfert 1 bv lOsterbrock fc Poggd 
(1985) on the basis that the FWHM of Ha line is less 
than 2000 km s~\ the ratio of [O III]A5007 to H 3 is 
less than 3, and they often show strong Fe II emission in 
4500-4680 A and 5105-5395 A regions. Though NLSls 
are identified by their optical properties they have even 
more remarkable properties in the X-ray band as com- 
pared to other Seyfert galaxies. They show evidence for 
strong excess of soft X-rays (dominant below ~ 2 keV) 
above the hard X-ray continuum extrapolation and rapid 
X-ray variability (doubling times of minutes to hours). 

Soft X-ray excess emission above a power-law contin- 
uum is usually identified as the steepening of the X- 
ray continuum below ~ 2 keV. Th is soft excess emis- 
sion was first o bserved by HE A O-l llPravdo et al ]ll981l) 
and EXOSA . T ijArnaiiri et~aTl IT985lTSmgh et a.llllflM . 
iBoller et all 1)19961) showed that AGN with steepest soft 
X-ray spectra in the ROS AT 'band tend to lie at the lower 
end of the H3 line width distribution and hence AGNs 
with soft excess emission are predominantly NLS1 galax- 
ies. Detailed X-ray study of NLSls was carried out with 
ASCA which confirmed the soft X-ray emission and, in 
addition, showed that the 2 — 10 keV continuum slope too 
is ste eper and anti-correlated with the FWHM of the Ha 
line iBrandt et al J 119971 iLeighlvl Il999b1: iVaughan et all 
1999). ASCA observations also showed that many NLSls 
have Fe K Q line arising from highly ionized iron, and 
hence the accret i on disks of NLSls must be ionized 
I pewangarJ 120021 iPorauet et.al]l2004t iBallantvne et all 
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12001 . NLSls are the most extreme X-ray variable ob- 
jects among the radio-quiet AGNs. NLSls frequently 
exhibit rapid (doubling time scale of a few hundred sec- 
onds) and/or l arge (up to a factor of 100) amplitude X- 
ray v ariability i)Boller et al]ll99ftlDewangan et alJl20Qll 
2002). The excess variance for NLSls is typically an or- 
der of magnitude higher than that observed for Broad 
Line Seyfert l's (BLS1). 

There are high-luminosity a nalogs of this class, the 
prototype of which is I Zw 1 (jPhillipsI fl976) . I Zw 1 
also has weak forbidden lines, narrow permitted lines, 
and strong Fe II lines. One of the most extreme nar- 
row line Fe II QSO is PHL 1092 (B=16.7 z=0 396) 
l|Bergeron fc KuntUll980l 11984 iKwan et. al.l 11995^ . Its 
Fe IIA4570/ii/ ( 3 ratio is 5.3, and i ts Hg line width is 
1800 km s" 1 <|Bergeron fc Kunthl Il98l . PHL 1092 
was rapidly variable during a ROS AT observation, and 
its soft X-ray spectrum was extremely steep, mod- 
eled by a power-law with photon index m ore than 4 
(|Forster fc Halpernlll996T: iBrandt et alJ ll999L From an 
ASCA observation lLeighlvll(T999b() found a huge soft ex- 
cess (which can be described as a black body) in the 
PHL 1092 spectrum which is 80% of the total flux. An 
evidence for a line at around 7 keV has also been found 
and this line could not be explained as an ionized iron 
line. 

I i < > • ' v iDewangan. fc Griffiths! (|2004fl made a de- 
tailed study of the XMM-Newton observations of I Zw 1 
and they showed that during the X-ray variability the 
power-law component varied without any change in 
the soft X-ray excess component. In contrast to this, 
PHL 1092 shows large variability and it has more than 
80% soft X-ray excess component. To understand the 
X-ray spectral variability of this high luminosity object 
we have analyzed an XMM-Newton observation on this 
source. We present results based on this observation, in- 
cluding flux selected spectroscopy of high and low flux 
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Fig. 1.— PN light curves of PHL 1092 in the 0.3 - 0.5 keV (top 
panel), 0.5 — 1 keV (middle panel) and 1 — 10 keV (bottom panel) 
bands with a bin size of 512 s. The high and low flux states used 
for flux selected spectroscopy are shown by dotted vertical lines in 
the top panel of the figure. 



states. 

2. OBSERVATION AND DATA REDUCTION 

PHL 1092 was observed by XMM-Newton observatory 
on 2003 January 18 using the European Photon Imaging 
Camera (EPIC) and the reflection grating spectrometer 
(RGS). The raw events were processed and filtered us- 
ing the most recent updated calibration database and 
analysis software (SAS v6 . . 0) available in 2004 March. 
Events in the bad pixels and those adjacent pixels were 
discarded. Only events with pattern — 4 (single and 
double) for the PN and 0-12 for MOS were selected. 
Examination of 16 s bin light curve extracted from the 
full field, excluding the source region, and above 12 keV 
showed high particle background in some observation in- 
tervals. The background flaring region is excluded using 
the criterion 'rate < 3' and minimum good time interval 
> 512 s for PN data. This resulted in 'good' exposure 
time of ~ 18.9 ks for the PN. And during these good time 
intervals the MOS1 and MOS2 background were reason- 
ably low, so we took the same time interval for MOS data 
also. The net count rate is 0.53 counts s _1 in PN and 
0.12 counts s- 1 in MOS1 and MOS2. 

3. TEMPORAL ANALYSIS 

X-ray light curves of PHL 1092 were extracted from 
the PN and MOS data using a circular region of 35" and 
45" respectively centered at the source position and in 
the two soft bands (0.3 — 0.5 keV and 0.5 — 1 keV) as 
well as in the hard band (1 — 10 keV). The energy bands 
were chosen to separate approximately the spectral com- 
ponents - a power-law (dominant at high energy) and 
a soft excess component gen erally observed from NLS1 
galaxies (see iLeiehlvl Il999bj) . Background light curves 
were extracted from source free regions with the same 
bin sizes, exposure requirements, as for the source light 
curves. The area of the region selected for the back- 
ground is four times the area of the region selected for 
the source. The background light curves were subtracted 
from the respective source light curves after appropriate 
scaling to compensate for the different areas of the ex- 
traction regions. 
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Fig. 2. — 2-10 keV PN spectrum and best fitted model of 
absorbed and red-shifted power-law. 



Figure 1 shows the PN light curves of PHL 1092 in the 
different energy bands with a bin size of 512 s. It is evi- 
dent that X-ray emission from PHL 1092 varied strongly 
during the XMM-Newton observation particularly in the 
soft bands (0.3 - 0.5 keV and 0.5 - 1 keV). There is a 
sharp decrease in the observed count rates starting from 
10 4 s (from start of the observation) and reaches a low 
level in about 5000 s. Two regions are marked as "high" 
(from start to 10 4 s) and "low" (1.5 x 10 4 s to 2.1 x 10 4 
s) in Figure 1 and demarcated by vertical dotted lines. 
A constant model fit to the complete data gives an un- 
acceptable value for the reduced x 2 (~ 7.7 and ~ 3.7 re- 
spectively in the two low energy bands). Fitting constant 
models to the data in two regions (namely the high-flux 
state and the low-flux state) gives much lower value of 
reduced \ 2 ~ 2.2 and ~ 0.73 respectively in 0.3 — 0.5 keV 
band and ~ 2.9 and ~ 0.87 respectively in the 0.5 — 1 keV 
band. The average count rate throughout the observa- 
tion in the 0.3 — 0.5 keV band is 0.28±0.01 s~ x , whereas 
in the high and low states the average count rates are 
0.36 ± 0.01 s _1 and 0.22 ± 0.01 s -1 , respectively. In the 
0.5—1 keV band the average count rates in the high 
and low states are 0.23 ± 0.01 s ~ x and 0.16 ± 0.01 s -1 . 
During the whole observation the average count rate was 
0.19 ±0.01 s _1 in this energy band. This shows that the 
observed count rates varied in a time scale of about 5000 
s at a very high significance level. Similarly the average 
count rate in the 1 - 10 keV band is 0.053 ± 0.005 s" 1 
in high flux state and 0.036 ± 0.006 s _1 in the low flux 
state, whereas the average count rate in this band during 
the whole observation is 0.047 ± 0.003 s _1 . So there is a 
clear variation of count rate in this band too. 

4. TIME-AVERAGED SPECTRAL ANALYSIS 

Photon energy spectra of PHL 1092 and associated 
background spectra were accumulated from the EPIC PN 
and MOS data using the source and background extrac- 
tion method described above. We have also analyzed the 
RGS data, but the observed count rates are too low to 
draw any definite conclusion. The source spectra were 
grouped such that each bin contained at least 15 counts. 
The EPIC responses were generated using the SAS tasks 
rmf gen and arf gen. All the spectral fits were performed 
with the XSPEC package (version 11.2.0: IArnaud1ll996j) 
and using the \ 2 -statistic. Unless otherwise specified, 
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Fig. 3.— The 0.3 -10 keV EPIC PN spectrum. The extrapolated 
power-law model is also shown. A strong soft excess is present 
below ~ 2 keV. 
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Fig. 4. — PN spectral data and the best-fit model consisting of 
two absorbed and red-shifted blackbody and a power-law (r ~ 2.1). 



the quoted errors on the best-fit model parameters were 
calculated for the 90% confidence level for one interesting 
parameter i.e., A% 2 = 2.7. 

4.1. 2 - 10 keV spectrum 

We began studying the spectrum by fitting a red- 
shifte d power-law model w ith a Galactic Njj of 3.53 x 
10 20 i|Murphv et. alJll996|) to the EPIC spectrum in the 
energy range 2-10 keV. This provides an acceptable 
fit giving a photon index T ~ 1.88 and a minimum 
X 2 (dof) — 22(27). The fit to the spectrum and the resid- 
uals are shown in Figure 2. There is no clear indication 
of any line feature near 6 keV. The non-detection of an 
Iron line is mainly due to the low level of detection sig- 
nificance above 5 keV. The upper limit for the equivalent 
width of a narrow iron line at 6.4 keV is 289 eV. 

4.2. 0.3-10 keV spectrum 

Extrapolation of the best fitting 2-10 keV power law 
to 0.3 keV shows a huge soft excess in the spectrum (Fig- 
ure 3). To fit the broad band spectrum, we used one red- 
shifted blackbody component to parameterize the soft 
excess. Adding a single blackbody provides a good fit 
X 2 ~ 192(196) to the soft emission but the power law 
slope steepens to 2.15. The blackbody temperature is 
found to be ~ 124 eV. Adding another blackbody gives 
an improvement with A% 2 ~ 7 for two extra parameters 
and the F-test null probability for adding this compo- 
nent is 0.027. The temperatures of the blackbody com- 
ponents are found to be 129 eV and 53 eV. The data and 
the folded model are shown in Figure 4. We have also 
tried using two red-shifted power-law which gives high 
reduced x 2 and is unacceptable. 

4.3. PN and MOS data 

We have made a joint fit to the PN and MOS data 
and find that the above model fits the data satisfacto- 
rily. The fit and the ratio of data to model are given in 
Figure 5 and the derived spectral parameters are given 
in Table 1. The ratio of data to model for the differ- 
ent instruments have been shown separately. In MOS1 
there is a slight deficit of counts at around ~ 1.4 keV, 
whereas in the PN this feature is less significant and in 
MOS2 this is not present at all. Examining the data to 



model ratio one can argue that there is a slight excess 
emission at high energies. But we did not get any im- 
provement in the fit by including line or edge parameters 
at the expected Fe lines and edges. To understand the 
nature of high energy emission properly we have put sev- 
eral markers (a,b,c,d,e,f) in Figure 5. Markers 'a' and 'b' 
are at 7.1 keV and 9.2 keV in the rest frame, respectively 
(corresponding to neutral and ionized edges of Fe). By 
examining the data to model ratio at these two energies 
one can conclude that there is no indication of edge in 
the data. Markers 'c' and 'd' correspond to 6.4 keV and 
6.9 keV (corresponding to neutral and H-like Fe emis- 
sion lines), and there is no apparent excess emission at 
these energies. The apparent line like emission feature 
in MOS2 data corresponds to unrealistically high energy 
(marker 'e' at energy ~ 8 keV in the rest frame). Again 
if we consider the edge like feature at high energy in the 
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Fig. 5.— PN. MOS1 and MOS2 spectral data and the best-fit 
model consisting of two absorbed and red-shifted blackbody and 
power-law (r ~ 2.1) (top panel). Ratio of data to model for PN, 
MOS1, and MOS2 are also shown separately for clarity. The labels 
a,b,c,d,e, and f correspond to the source rest frame energies of 7.1 
keV, 9.2 keV, 6.4 keV, 6.9 keV, 8 keV and 10.5 keV, respectively 
(see text). 
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TABLE 1 

Best-fit spectral parameters of PHL 1092. 



Model 



model 1 



Component 



Parameter 



PN 



PN+MOS 



model 2 
PN 



Blackbody 
Blackbody 
Power law 
Comptt 

Comptt 



Total 



kT (eV) 
n-bb 2 

f bb (0.3-10 keV) 3 
kT (eV) 
n-bb 3 

/ w (0.3-10 keV) 4 

r 

/ pi (0.3-10 keV) 3 
fcT seed (eV) 

kTpi asma (keV) 
T 

n c 5 

/ c (0.3 - 10 keV) 3 
kT seed (eV) 



f c (0.3 - 10 keV) 3 

f obs (0.3-10 keV) 3 
f intr (0.3-10 keV) 3 
L intT (0.3-10 keV) 6 



129.0^1^ 
6-38±°-.|| 
9.91 

p-o + 13.3 
0^-«_24.5 
-.o +217.6 
ro.U_ 11 2 

1.12 

1+0.21 
-0.16 



130.1] 

6.04 4 



2.04 1 



1.19 



+0.31 



0.19 

2.95 



+3.4 
-2.8 
,+0.40 
-0.60 

9.51 
54.3± 3 ; 2 . 
13.8± 37 2 ° 

1.39 

9 ic+0.13 
z - lo -0.15 
-, o 7 +0.18 
L - 3 < -0.19 

3.06 



110 
4 

0.10 
1.4 



+10 
12 
+ 12.0 
-4.1 
+0.39 
0.09 
+ 104.2 
0.9 
11.40 



0.32 
59.2 

0.51 
0.002 
2.26 



+0.25 
-0.32 
16.-1.3 
57.2 
+5.05 

+B°870 
0.001 



X 2 



14.0 
20.8 
20.1 

186/194 



14.0 13.6 

20.9 19.9 

20.2 19.2 

334/326 189/191 



1 Parameter values for spectral fits EPIC PN and MOS data. 
Spectral models were modified by the Galactic absorption. 



3.53 x 10 2U cm" 2 ), model 1: WABS(BB+BB+PO), 



model 2: WABS(COMPTT +COMPTT) where WABS: photoelec- 
tric absorption, BB: blackbody PO: power-law, COMPTT: Comp- 
tonization 

2 Blackbody normalization in units of 10 — 5 X 10 39 erg s — 1 /(dio) 2 , 
where d is the distance in the unit of 10 kpc 
3 Flux in the unit of 10 -13 erg cm -2 s _1 

4 Powcrlaw normalization in units of 10 -4 photons cm -2 s _1 
keV" 1 at 1 keV 

5 Comptonization normalization in units of 10 _3 photons cm -2 
s- 1 keV" 1 

6 Source luminosity in the unit of 10 44 erg s _1 



data the corresponding energy will be ~ 10.5 keV 
(marker £ f in the figure) in the rest frame which is unre- 
alistically high. Hence we can conclude that it is not pos- 
sible to quantify the high energy excess by these known 
physical phenomenon. One possible reason of this ex- 
cess emission at around 8 keV or deficit of counts at 
around 10.5 keV is an artifact of incorrect background 
subtraction. In this range background is dominated over 
the source counts giving rise to an apparent curvature. 
But if the curvature is true and due to some physical 
phenomenon then line emission from relativistic ionized 
accretion disk (laor) cannot be ruled out. By putting a 
laor line we get an improvement of A% 2 ~ 8 (F-test null 
probability ~ 0.2), and the line energy is ~ 6.7 keV in 
the rest frame (see section 6.4) for five extra parameters. 

4.4. Two component thermal Comptonization model 

The red-shifted and absorbed blackbodies and power- 
law model described the data adequately, but the tem- 
perature of the blackbodies (kT ~ 130 eV) and (kT ~ 50 



eV) are unrealistically high for a standard accretion disk 
(see section 6.3). A corona with two electron popula- 
tions with distinct temperatures can give rise to both 
a soft X-ray excess and a power-law. We investigated 
if such a process is a viable model for the X-ray spec- 
trum of PHL 1092. The contribution of this soft excess 
component to the 0.3 — 2 keV band unabsorbed emis- 
sion is ~ 80%. As the temperature of the soft excess 
component is much higher than that expected from an 
accretion disk, this component is unlikely to be the in- 
trinsic disk emission. Replacing the soft component by 
comptt we get a good fit (x 2 — 190 for 193 dof). Then we 
replace the power-law with a hotter comptonized compo- 
nent. The fit is as good as previous (\ 2 = 189 for 193 
dof), since the EPIC instrument cannot determine the 
value of the electron temperature if kT > 4 keV. The 
exponential roll-over of the comptonized component oc- 
curs at ~ 4&T, and below this point the emission appears 
as a power-law. Thus if AkT is very much greater than 
■~ 10 keV, the temperature cannot be constrained over 
the XMM-Newton band, and the fit is indistinguishable 
from a power-law over the same energy range. Because 
of this the resultant % 2 value for the high temperature 
Comptonization fit is virtually identical with the power- 
law fit. The parameter values of two comptonization fit 
are given in Table 1. 

5. FLUX-SELECTED SPECTROSCOPY 

To investigate whether the flux variability is due to 
changes of the spectral parameters we have carried out 
a spectral analysis at low and high flux levels (shown in 
Figure 1 by dotted vertical lines). We extracted two av- 
erage spectra covering - 10 ks (high flux state) and 15 - 
21 ks (low flux state) in the 0.3 - 10 keV band. The aver- 
age PN count rates are 0.64 ± 0.02 counts s _1 (high flux 
state) and 0.41 ± 0.02 counts s _1 (low flux state). We 
have used the model (absorbed and red-shifted power- 
law, and two blackbody) for both the spectra. First we 
fixed all the parameter previously found values for the 
combined fit of PN and MOS data and obtained a value 
of x 2 (dof) of 242(141) and 285(85) for the high and low 
flux data, respectively. These values improved to 165 
and 77 when the relative normalization of the two data 
sets are changed. The fit improved marginally (x 2 (dof) 
of 152(144) and 76(87), respectively) when all the nor- 
malizations are allowed to vary mutually. The observed 
count rates are too low to investigate spectral variability 
at a greater detail. In Table 2 the fit parameters and 
the observed fluxes are given. It can be seen from the 
table that the source is highly variable during the obser- 
vation. The change in luminosity is ~ 35% between the 
two states. 

6. DISCUSSION 

This paper presents an analysis of ~ 20 ks XMM- 
Newton observation of the narrow- line QSO PHL 1092. 
The time-averaged X-ray spectrum of PHL 1092 shows 
a power-law continuum, a soft X-ray excess described by 
two red-shifted blackbody components below ~ 2 keV. 
Above 2.0 keV, the power-law component has a slope 
of r ~ 2.1 which is steeper than that of most BLS1 
galaxies (r ~ 1.9) but is similar to other NLS1 galax- 
ies e.g, Ton S180 (r = 2 26±£ f 7 - IVaughan et alJ f2002L 
Mrk 335 (r = 2.29±0.02: lGondoin et aDl2002]) . Akn 564 
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TABLE 2 

Best-fit spectral model parameters derived from the 
spectral fitting of the pn spectra of phl 1092. 







Flux States 1 


Component 


Parameter 


(high) 


(low) 



BB (kT- 
BB (kT- 
PO 


-130 eV) 
-50 eV) 


n b b 2 
n bb 2 


6.87l°;»? 

24 6+ 6 7 
-, Qo+0.13 
1 '°°-0.13 


4 59+ 43 
74+° ;74 


Total 




f obs (0.3-10 keV) 4 
f intr (0.3-10 keV) 4 
L inU (0.3-10 keV) 5 
X 2 min /dof 


16.01 

24.53 
23.64 
152/144 


10.07 
14.91 
14.37 

75/88 



erg s 



V(dio) 2 



High: High flux state; Low: Low flux state 
2 Blackbody normalization in units of 10 — 5 X 10 39 

where d is the in 10 kpc distance. 
3 Power-law normalization in units of 10 — 4 photons cm" 2 

keV" 1 at 1 keV. 
4 Flux in units of 10 — 13 erg cm" 2 s" 1 
5 Source luminosity in units of 10 44 erg s — 1 



(r ~ 2.50 - 2.55: lVignali et alJl2003|) . The X-ray contin- 
uum steepens at lower energies, resembling the soft X-ray 
excess emission generally observed from NLS1 galaxies. 
The strength of this soft X-ray emission is quite high in 
PHL 1092 compared to other NLSls. The contribution 
of the soft X-ray excess emission to the 0.3 — 2 keV band 
X-ray emission is about 80% in PHL 1092. Such high 
soft excess has been reported earlier bv ILeighlvl (^999b) 
based on the ASCA data. Variability at similar variabil- 
ity tim e scale has been reported by iForster fc Haloernl 
1)1996(1 based on the ROSAT data. 

6.1. Black hole mass for PHL 1092 

Before we discuss the results in detail, we estimate 
the mass of the black hole for PHL 1092. Calcula- 
tions of the black hole mass rely on the assumption 
that the dynamics of the optical broad-line region gas is 
dominated by the gravitational potential of the central 
super-massive black hole. To calculate the black hole 
mass for PHL 1092, we used the results of Kaspi et al. 
(2000), based on a reverberation study of 17 quasars. 
iKaspi et alJ ((2.000) have determined an empirical rela- 
tionship between the size of the broad-line region (Rblr) 
and the monochromatic luminosity (XL\) as follows. 



Rblr = (32.9; 



-2. On, 
-1.9J 



AL A (5100) 



10 44 erg s" 



0.700±0.033 



It — days 
(1) 

We calculated the flux approximately from the IDS spec- 
trum of PHL 1092 obtained from iBergeron fc Kuntbl 
(1980) at the rest wavelength of 5100 A (A/a (5100 A) 
~ 3 x 10 -12 erg cm" 2 s" 1 ). The mass of the black hole 
is given by Mbh — rv 2 G~ x , To determine v, the ve- 
locity, we correct Vfwhm °f the Ha emission line by a 
factor of v3/2 to account for velocities in three dimen- 
sions . We us e d vfw hm = 1790 km s _1 for the H3 line 
from ILeighlvl l)1999a|) . The mass is then, 



M = 1.464 x 10 6 



RBLR \ ( V FWHM \ 2 { ) 

It -days ) Vl03kms-iJ ° 1 j 

Using the cosmology go = and Hq — 50 km 
s~ 1 Mpc~ 1 , this gives us a value f or the central black 
hole mass of - (1.6 ±0.6) x 10 8 M^. ICzernv et al.l<|2001f) 
have reported that the mass of black hole in PHL 1092 
is 1.2 x 10 6 M Q from X-ray variability studies and it is 
1.8 x 10 8 Af Q from accretion disk fitting. 

6.2. Rapid variability and the efficiency limit 

PHL 1092 showed rapid and large amplitude X-ray 
variability during the XMM-Newton observation. X-ray 
emission from PHL 1092 decreased by ~ 35% in 5000 
s, corresponding to a change in the 0.3 — 10 keV lu- 
minosity of AL ~ 10 45 erg s" 1 . This is a remarkable 
variability for a radio-quiet quasar as the time scale is 
a factor of — 10 shorter than the dynamical time scale 
for a 10 s Mq Schwarzschild black hole but comparable 
to the dynamical time scale of Kerr black hole of same 
mass. Such variability is rare among radio-quiet quasars. 
Other radio-quiet quasars to sho w such rapid and large 
ampl itude variability are I Zw 1 ij Dewangan, Griffiths 
2004), RX J1334.2+3759 JDewangan et al.ll2002t) . 

The major consequence of rapid and large amplitude 
variability is that the source has to emit very efficiently. 
It is possible to calculate the efficiency of conversion of 
rest m ass into e nergy using simple arguments originally 
due to iFabianl l(l979T) . If the X -ray emitting region is 
an accreting spherical cloud of ionized hydrogen, then 
there is a limit on the radiative efficiency, rj, given by 
77 > 4.8 x 10 -43 AL/Ai, where AL is the change in lu- 
minosity in the rest-frame time interval At. The change 
in luminosity between the two states is ~ 5 x 10 44 ergs 
s" 1 in 5000 s. Assuming uniform and spherical emis- 
sion this gives a lower limit to the efficiency, rj > 0.1. 
This efficiency can be compared with the maximum effi- 
ciency possible for a Schwa rzschild black hole (77 ~ 0.06; 
IShapiro &: Teuko lskvllT9 83) and a Kerr black hole (77 ~ 
0.3: lThorndll974]) . Thus the rapid variability observed 
from PHL 1092 is consistent with that expected for a 
Kerr black hole. 

6.3. The soft X-ray excess emission 

Although the soft X-ray excess emission above a power- 
law continuum is well described by a blackbody compo- 
nent, its temperature is too high for the intrinsic emission 
from an accretion disk. There is a relationship between 
the disk temperature of a standard thin accretion disk, 
accretion rate a nd the mass of the black hole, given by 
l|Petersonlll99H) . 



T(r) 



6.3 x lOW^Mg" 174 f — 



-3/4 



K 



(3) 



where r$ is the Schwarzschild radius, and Mg is the 
black hole mass in units of 1O 8 M0. The temperature 
at the last stable orbit (r = 3rs) for a black hole is 
— 22 eV for a black hole mass of 1.6 x 1O 8 M0 accret- 
ing at the Eddington rate. This is an upper limit to the 
disk temperature, because m < 1, and outer regions of 
the disk will have lower temperatures. This can be rec- 
onciled by assuming that the low energy excess is due 
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Fig. 6. — Unfolded PN spectral data and best fit model com- 
ponents: two absorbed and red-shifted blackbody and power-law 
(above) and two absorbed and red-shifted blackbody, power-law 
and a laor line (below). 




0.5 1 2 5 

Energy (keV) 

Fig. 7.— PN, MOS1, MOS2, ASCA-SISO, ASCA-SIS1 spectral 
data and the best-fit model consisting of two absorbed and red- 
shifted blackbody and power-law (top panel). ASCA-SIS data to 
model ratio freezing the power-law index of ASCA data to that 
found in XMM-Newton data (middle panel) and the same model 
but keeping the value of power-law index of ASCA data free (lower 
panel). PN, MOS1 and MOS2 data to model ratio have not been 
shown for clarity. 



to a population of electrons very close to the black hole 
giving a Comptonized spectrum. In fact the observed X- 
ray spectrum in the 0.3 - 10 keV region can be fit with 
a two component Comptonization model. An alterna- 
tive viable solution is the 'slim disk' l)Abramowicz et aD 
119881: iChen fc'W ang 2003, an d references therein) which 
allows temperature higher than the standard disk, and 
takes place at high accretion rates. 

6.4. Possible evidence for a high-energy curvature and 
absorption 

Recently iGallo et. al.1 l)2004|) reported a high energy 
curvature around ~ 7 keV in the same data of PHL 1092. 
They interpreted the high energy curvature as due to 
light bending effects very close to the black hole. When 
we analyze the data using the latest analysis tools (SAS 
6.0.0) and latest calibration, we find very little evidence 
for a high energy curvature and the curvature is cen- 
tered around 7.5 keV. A narrow Gaussian line cannot 
give satisfactory fit to the high energy curvature. How- 
ever a broad Gaussian (width ~ 1 keV) can give an visual 
improvement of the spectrum (but not statistically sig- 
nificant, A% 2 ~ 3 for three extra parameters, F-test null 
probability ~ 0.4) but the line energy becomes greater 
than 7.2 keV which is clearly unphysical. However a line 
from a disk around a rotating black hole (Laor 1990) can 
be used to handle the high energy curvature. The fit is 
marginally better than the previous fit (see section 4.3). 
The unfolded spectrum without the Laor line is given in 
the top panel of Figure 6 and with the line is given in 
the bottom panel. The line energy is ~ 6.7 keV and the 
equivalent width is ~ 2.5 keV. Hence we can conclude 
that Gaussian lines are not required to fit the data, but 
Laor lines from ionized accretion disk cannot be ruled 
out. The large equivalent width of the line seems unphys- 
ical and in contrast with the fact that the line equivalent 
width diminishes with increasing continuum luminosity 



(Baldwin effect). IGallo et. all l)2004|) also reported an 
absorption line around ~ 1.4 keV. This absorption fea- 
ture was not present during the 1997 ASCA observation 
(Leighly 1999b). In section 4.3 we discussed why inclu- 
sion of an absorption line or an edge is not feasible. The 
major reason for the difference in the results given in 
Gallo et al. (2004) and the present work is the use of 
most recent calibration and SAS version. 

To further understand the high energy spectral shape, 
we h ave reanalyzed t he ASCA data observed on 1997 
,Tulv iLeighlvll999albl) . The ASCA-SISO and ASCA-SIS1 
data are simultaneously fitted with the XMM-Newton 
EPIC-PN and EPIC-MOS data. To avoid the calibra- 
tion uncertainties of ASCA-SIS (see ASCA calibration 
page) below 1 keV we take 1-10 keV energy band. We 
fixed all the parameters obtained from the joint fit of 
EPIC-PN and EPIC-MOS. The relative normalizations 
between the instruments are allowed to vary. The rel- 
ative normalizations of ASCA-SISO and ASCA-SIS1 are 
1.51 and 1.38 respectively. The fit obtained is quite sat- 
isfactory (x 2 ~ 397 for 377 dof), indicating that over 
the different epochs of observation variability can be ex- 
plained by overall flux variation rather than spectral pa- 
rameter variation. The observed and fitted spectra are 
shown in Figure 7, along with residuals (plotted as ra- 
tio of the data to the model in the bottom panel of the 
figure). It can be seen that there is a reasonable agree- 
ment between the data and the model. Changing the 
power-law index improved the fit with a A\ 2 of 17, with 
the power-law index changing to 1.73. In this case the 
high energy excess is even lower (see Figure 7 bottom 
panel). We can conclude that once the low-energy spec- 
tral shape is appropriately modeled, there is no necessity 
of invoking any high energy curvature. The data qual- 
ity, however, is insufficient (particularly above 5 keV) 
to make any definitive statements about other possible 



XMM- Newton observation of PHL 1092 
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spectral features like edges, reflection features etc. 

6.5. Comparison with the high-soft state of Galactic 
black hole sources 

All the available data suggest a smooth spectrum with 
all the spectral features varying at a fast time scale. The 
time scale of variability, however, corresponds to a physi- 
cal size (light crossing time) of 1.5 x 10 9 km, which corre- 
sponds to only 3 Schwarzschild radius for a 1.6 x 10 8 M© 
black hole. Since both the spectral components appear 
to be varying in this time scale, we can conclude that 
the X-ray emission arises from very close to a rotating 
black hole. Such emission mechanism is seen in the high- 
soft states of Galactic black hole (GBH) binary sources. 
The accretion disk reaches very close to the black hole 
and bulk motion near to the b lack hole can give rise to 
the o bserved rapid variability l)McClintock fc Rem illard 
12004(1 . The rapid variability observed in PHL 1092 on 
time scale < 10000 s can be com pared with the millisec- 
ond f la res observed in GBHs (e .g. TZdziarski fc Gierlinskil 
12004(1 . iMcHardv et. all l|2004fi have compared the vari- 
ability of NGC 4051 in long and short time scales with 
that of Cygnus X-l and have concluded that the X-ray 
variability in NLS1 galaxies scale better with Cyg X-l in 
its high state. 

The observed 0.3 - 10 keV X-ray luminosity corre- 
sponds to 7% (or 10% after correcting for Galactic ab- 
sorption) of the Eddington luminosity (LEdd)- If X- 
ray luminosity represents a fraction (typically ~ 10%) 
of bolometric luminosity then the bolometric luminos- 
ity of PHL 1092 will be ~ L Edd - The high-soft states of 
Galactic black hole sources show a total X-ray luminosity 
above ~ 0.1 L Edd l|Zdziarski &; GierlinskilfcOO^ . 

For the fast correlated X-ray variability, we can have 
the following physical picture. The low energy black- 
body arises from a optically thick disk extending up to 
the last stable orbit of a rotating black hole. Within 
this radius, matter gets accreted at a supersonic velocity 
which can produce bulk-motion Comptonisation mani- 
fested as the power-law. Non-thermal activity (by the 
possible magnetic field generated by the equipartition of 
energy) can generate high energy electrons which would 
be cooled by Comptonization by the copious amount of 
low energy photons from the low-energy blackbody emis- 
sion. This can produce highly saturated Comptonization 
spectrum modeled as a blackbody in the present work. 
We invoke the bulk-motion Compto nization model given 
bv iChakrabarti fc Titarchukl ll 19951) to estimate the rel- 
evant energetics for such a simplified model. The to- 
tal luminosity of the ~ 50 eV blackbody corresponds 
to ~ 5.4 x 10 45 ergs s _1 which is the expected lumi- 
nosity from a thick accretion disk extending up to 3 
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Schwartzchild radius for a 1.6 x 10 8 M rotating black 
hole with 0.2 Eddington accretion rate. The power-law 
luminosity is about 10% of the total luminosity, which is 
similar to t hat seen in Galactic black hol e sources in their 
high state (Zdziarski & Gicrlinskil l2004D . Comparison of 
the timing and spectral properties of such high luminous 
NLS1 galaxies with the high soft states of Galactic black 
hole sources would provide additional handle to constrain 
such a model. 

7. CONCLUSIONS 

We presented spectral characteristics of the narrow line 
QSO PHL 1092 based on an XMM-Newton observation. 
The main results are as follows. 

1. The 0.3 - 10 keV spectrum of PHL 1092, ob- 
tained with XMM-Newton, consists of three intrin- 
sic spectral components, namely a steep (Fx ~ 2.1) 
primary continuum described by a power-law, a 
soft X-ray excess component described by two red- 
shifted blackbody (kT ~ 130 eV and kT ~ 50 eV). 
Iron fluorescent lines are not detected, primarily 
due to the low detection level at > 6 keV. The in- 
trinsic photoelectric absorption component is con- 
sistent with the Galactic absorption. 

2. PHL 1092 showed a significant X-ray variability 
with changes in the luminosity AL ~ 10 45 erg s _1 
on a time scale of ~ 5000 s. The radiative efficiency, 
7] > 0.1, inferred from the variability is consistent 
with X-ray emission outside the last stable orbit 
around a Kerr black hole. 

3. The long (comparing ASCA and XMM-Newton ob- 
servation) and short term variability of PHL 1092 
can be explained by change in overall normalisation 
without any spectral changes. 

4. The existing data do not require any high energy 
curvature, but we cannot completely rule out the 
existence of a red-shifted line from the accretion 
disk near a fast rotating black hole (Laor line). 
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